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We calculate the effect of Coulomb interactions on the 
average density of states vluj) of two-dimensional disordered 
electrons. It is shown that for weak disorder the most singular 
terms in the perturbative expansion of v{u)) can be summed 
by means of a simple gauge transformation, which also es- 
tablishes a relation between the low-frequency behavior of 
v{lo) and the average conductivity oiuj). Using this relation, 
we show that if livnuj^o o-{u>) < oo, then u{u!) ^ C|t<;|/e* for 
— > 0, where C is a dimensionless constant and e is the 
charge of the electron. This implies that a normal metallic 
state of disordered electrons in two dimensions is not a Fermi 
liquid. 

PACS numbers: 71.23.-k, 71.30.-Fh, 72.15.Rn 



In two dimensions arbitrarily weak disorder causes 
non-interacting electrons to localize, so that the conduc- 
tivity vanishes in the thermodynamic limit 0| . For many 
years it was generally accepted that electron-electron in- 
teractions do not qualitatively change this scenario. The 
recent experimental discovery of a metal-insulator tran- 
sition in two-dimensional {2d) semiconductor devices |^] 
was therefore a surprise, and showed that the current the- 
oretical understanding of electron-electron interactions in 
disordered electronic systems is incomplete. The experi- 
ment Hwas analyzed in the light of a generalized scaling 
theory p], but a microscopic understanding is still lack- 
ing. Motivated by these exciting new developments, in 
this work we shall re-examine the effect of long-range 
Coulomb interactions on the average density of states 
(DOS) I'iiij) of disordered electrons in two dimensions. 

It is well known that Coulomb interactions drastically 
modify the DOS in the vicinity of the Fermi energy fi. 
For example, in the strongly localized regime they give 
rise to a gap in the DOS at the Fermi energy (Coulomb- 
gap), which in 2d is of the form W 



ro = [(27r)2l?oz.o] 



Tl = ToIVqk'^TqY 



(3) 



zy(ijj) ~ C\uj\/e^ , Lu 
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where C is a dimensionless coefficient and e is the charge 
of the electron On the other hand, for weak disorder 
(i.e. for kp£ ^ 1, where kp is the Fermi wave- vector and 
£ is the elastic mean free path) a perturbative calcula- 
tion to first order in the dynamically screened Coulomb 
interaction yields in two dimensions , 



where we have defined 



l-^ln(|c.|ri)ln(|.;|ro) 



(2) 



Here i/q — m/27r is average DOS (per spin) at the Fermi 
energy of non- interacting electrons with mass m, Vq = 
vp£/2 is the classical diffusion coefficient (where vp is 
the Fermi velocity), tq = £/vp is the elastic lifetime, and 
K is the Thomas-Fermi screening wave- vector Q]. Note 
that ro — {nkpi)^^ is a dimensionless measure for the 
resistance of the system at the frequency scale w ~ Tq^. 

Obviously the correction term in Eq.(|^) diverges for 
LU 0, indicating the breakdown of perturbation theory. 
It is tempting to speculate that at sufhciently small fre- 
quencies Eq.(||) should eventually cross over to Eq.(|]), 
at least if the localization length ^ is finite. Because for 
non-interacting electrons in the weakly localized regime 
^ (X exp[l/ro] is exponentially large, the crossover should 
only be visible at exponentially small frequency scales. 
Indeed, resummations of the leading logarithms |P-pd|] in 
the regime of weak disorder (where ro <C 1) indicate that 
in 2d the DOS vanishes for w — > 0, but the precise manner 
in which this happens has not been clarified |jl^,|l2| . 

It is instructive to recall the origin of the In^-divergence 
in Eq.(|^). The average DOS of an interacting Fermi sys- 
tem can be calculated from the Fourier transform of the 
retarded Green's function Gfl,(r,r',t), 
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iy{Lu) = Im / die'"*Gi?(r,r,t) , 

TI" J-oo 
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where the overline denotes averaging over the disorder. 
By translational invariance, the average is independent 
of the position r. In a perturbative approach, the cor- 
rections to the DOS are obtained by expanding Gr in 
powers of the interaction and averaging each term in 
the expansion over the disorder, using the impurity dia- 
gram technique. The Feynman diagram responsible for 
the In^-singularity in Eq.(^) is shown in Fig.|l](a). This 
diagram has the special property that at small wave- 
vectors (|q| <C £~^) and low frequencies (|a;| <C Tq^) the 
diffusive motion of the electrons leads to singular cor- 
rections to the charge vertices, which in turn strongly 
enhance the contribution from the long-wavelength and 
low-energy part of the dynamically screened Coulomb 
interaction. The vertex corrections involve so-called 
Diffuson-diagrams (Fig.||(b)), which within the Matsub- 
ara formalism renormalize the charge vertices for frequen- 
cies in the regime <I;„(a)„ -I- tUrn) < by a singular factor 



A(q, iwm;«w„) = [ro(I?oq^ 



(5) 



Here ujn — 27r(n-|-i)//3 is the fermionic frequency carried 
by the external Green's functions, and tOm = 2iTm/(3 is 
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the bosonic frequency carried by the interaction. (3 is the 
inverse temperature. Due to these singular vertex correc- 
tions, the main contribution to the diagram in Fig|^(a) 
comes from energy-momentum transfers satisfying 

\uJm\/VQK < |q| < {\uj„,\/Vo)^'^ , \uJra\ < T^^ ■ (6) 

In this regime, the dynamically screened averaged 
Coulomb interaction is within the random-phase approx- 
imation (RPA, see Figj^(c,d)) given by 



RPA 



(7) 



From Eqs.(^-|7|) it is now straightforward to show that 
the exchange diagram in Fig.^a) indeed gives rise to the 
In^-correction in Eq.(^). 

For our calculation described below it is crucial that 
this In^-correction is entirely due to Coulomb interactions 
with momentum transfers small compared with kp. In 
other words, the singularity is due to forward scatter- 
ing. Recently non-perturbative methods have been de- 
veloped to study electron-electron interactions with 
dominant forward scattering in clean systems. We shall 
now show that the functional bosonization approach de- 
scribed in Ref. ||l3|(a) can be generalized such that for 
weak disorder the dominant corrections to the DOS can 
be summed to all orders in the dynamically screened 
Coulomb interaction. 

We begin with a number of exact manipulations of the 
real space, imaginary time Green's function G(r, r',r — 
r') for a given realization of the disorder. Representing 
G as a Grassmannian functional integral, decoupling the 
interaction in the particle-hole channel by means of a 
Hubbard-Stratonovich (HS) field ^and finally integrat- 
ing over the fermions, we obtain [|l3[ (a) 



G(r,r',r 



e(r,r', r,r') 



/I?{0}e-SW} 



where the effective action is given by mA 



- Trln[l - GoV] 



(8) 



(9) 



and CJ(r, r', T, t') satisfies 

[-dr + Vl/2m + n- V{y, t) - U{v)\ ^(r, r', r, r') 

= 5{v-v')5*{t-t') . (10) 

Here V is the volume of the system, /q = 27re^/|q| is 
the Fourier transform of the bare Coulomb interaction in 
2d, and V and Go are infinite matrices in momentum- 
frequency space. The matrix elements of V arc [yjfefe' = 
il3~^(l)h-k' J and [Go]/cfc' is the non-interacting Matsubara 
Green's function for a given realization of the disorder. 
In Eq.(p^) U{r) is the (short-range) disorder potential, 
y(r,r) = ip-'E.K and 5*{r) = E„e*"""- 

Our strategy is now to separate the dangerous modes 
of V(r, t) with energy-momentum transfers in the regime 



(^ (which are responsible for the In^-singularity in 
Eq.(y)) from the harmless modes outside this regime. Let 
us therefore define V{r,T) — Vd(r,T) + Vh{r,T), where 
Vd(r, r) — ifi~^^q 4>q, and the prime indicates that the 
sum is restricted to the regime (||). We then substitute 
into Eq.(p^) the ansatz 

(?(r,r',r,r') =C?i(r,r',r,T')e*("-^^)-*("-''^') , (11) 

and seek a solution where ^(r, r) depends only on the 
dangerous modes Vd(r, r). It is not difficult to show that 
one obtains an exact solution of Eq.(pX|) by choosing Qi 
and $ such that -9r<I>(r, r) = Vd(r,T) and jl^ 

[-dr - i-iVr + A(r, r))V2m + /i - F,,(r, r) - ;7(r)] 
X ^1 (r, r', T, r') -<5(r-r')r (r-r') , (12) 

where A(r, r) = — iVr<I'(r, t) is a longitudinal vector po- 
tential. The equation for <I>(r,T) can be solved trivially, 
and we obtain for the gauge factor in Eq.(ni]) 



^^(rjr)— <i>(r',r') 



I /- y^' Vq 

I (3^ iUrn 



gi(q-r-w„T) _ gi(q-r'-ij„,T' 



(13) 



Nested exponentials of this type are familiar from 
bosonization p^(a). The crucial observation is now that 
Eq.(|l3|) contains the leading \t? -singularities to all or- 
ders in perturbation theory. To see this, note that in 



Eq.(|l2[) we have succeeded to eliminate from Eq.(lO) the 
dangerous part Vd(r, r) of the scalar potential in favour 
of a longitudinal vector potential. Although the associ- 
ated current vertices can still be dressed by Diffusons, 
the vertex corrections are less singular than in the case 
of density vertices [Q. For example, it is not difficult to 
check that the diagram analogous to Fig.|l|(a) with cur- 
rent vertices instead of density vertices does not give rise 
to any In^-singularities. From Eq.(pT|) it is also obvious 
that the gauge factor e*(''''^)^*('' ''^ ) cancels in gauge in- 
variant correlation functions such as the polarization or 
the conductivity, so that we immediately see that In^- 
singularities do not appear in the perturbative calcula- 
tion of these quantities ^,|| . 

To make further progress, we have to make sev- 
eral approximations, all of which can he systemati- 
cally improved. (a) First of all, because the fields 
A(r,r) and T4(r,r) in Eq.(p^) do not generate any In^- 
singularities in perturbation theory, let us approximate 
C/i(r, r', r, r') « Go(r, r', r — t'). This is sufficient to ob- 
tain the leading infrared behavior of the DOS, which is 
known to be dominated by In^-singularities g. Eq.(||) 
can then be written as G(r,r',T) = Go(r, r', r)e'5''"^'' '"^^ 
where Q(r,r',r) = ln(e*(""'^)-*(""''0)). Here (...) de- 
notes averaging over the HS-field as defined in Eq.(||). 
The Debye- Waller factor Q{t,v',t) can be calculated 
systematically via a linked cluster expansion in powers 
of the RPA interaction /q^^^^^^ for a given realization 
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of the disorder. For the average DOS we only need 
(5(0, 0,t) = Q{t). (b) We now calculate Q{t) in Gaus- 
sian approximation, retaining only the first order term in 
the linked cluster expansion. A simple calculation yields 
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E 

qq'w„ 



RPA 
qq','ia;„ 



■[1 - cos(cj,„r)] 



(14) 



In clean systems with dominant forward scattering the 
Gaussian approximation can be justified via Ward- 
identities |13|| . In the presence of disorder the situation 
is not so simple. We shall come back to this problem 
below, where we shall argue that, at least as long as 
the system has a finite static conductivity, Eq. ([l^) leads 
to qualitatively correct results for z^(a;). Note that by 
integrating over the HS-field before averaging over the 
disorder we have eliminated the disorder-dependent de- 
nominator in Eq. (H) . (c) Finally, we approximately aver- 
age over the disorder, using the factorization G(r, r, t) « 
G'o(r, r, r) exp[Q(r)]. Diagrammatically this amounts to 
ignoring all terms where impurity lines connect polariza- 
tion bubbles to other parts of a given Feynman diagram, 
a standard approximation in the usual diagrammatic ap- 
proach to weakly disordered electrons ||]. Using Eqs.(||) 
and (|^), we obtain in the limit V — > c» and f3 —^ oo 

W) = \ r —^0 Hu/V^k'') [1 - cos(c^r)] . (15) 



For large t this imphes Q{t) ^ — ^ ln(T/ri) ln(T/To), 
where ro and ti arc defined in Eq.(g). From Eq.(^) we 
finally obtain after analytic continuation 
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- sin(^ 
at exp 



-^ln(i/ri)ln(i/ro) 



(16) 

Expanding the right-hand side of Eq. ( [iq) to first order in 
ro, we recover the perturbative result (0), which is valid 
as long as — ln(|cij|ro) ^ \f\JrQ. Keeping in mind that 
1/ro 3> 1, we find that there exists an intermediate fre- 
quency range ■\/l/ro 1^ — ln(|u;|To) 1^ 1/ro where Eq.(|l^) 
can be approximated by 



v(ijj) K. i/Q exp 



-^ln(|w|ri)ln(|tj|ro) 



(17) 



in agreement with Finkelstein |Q. This is a non-trivial 
check that the gauge factor (|l3|) indeed contains the dom- 
inant singularities. Note that we have derived Eq.(|l7|) 
without using replicas. For frequencies below the expo- 
nentially small scale exp[— 1/ro] Eq.(p^ is not valid 
It is easy to show, however, that in this regime 
Eq.(|16[) smoothly crosses over to the linear Coulomb 
gap given in Eq.([I]), with the numerical constant C = 
47r~^/^Y^exp[l/ro]. The amazing fact is that for a; — > 
the prefactor vq in Eq. ( p^ ) disappears and is replaced by 
\uj\/e*, which in 2c? has the same units as I'o. 



We now argue that if the conductivity a{uj) does not di- 
verge for w ^ 0, then the true asymptotic low-frequency 
behavior of the DOS is indeed given by Eq.(|l|), although 
our above result for the numerical value of C is not re- 
liable. Let us therefore recall that the parameter ro in 
Eq. (|l5|) is the dimensionless resistance at frequency scale 



Tq . With the above approximations (b) and (c) we have 
effectively replaced S{(t)} in Eq.(||) by its Gaussian aver- 
age, which in the regime (0) can be written as 



V 



ro\uj„ 



(18) 



In deriving Eq. (|15|) , we have assumed that ro is not renor- 
malized for w — > 0, which is in general incorrect: higher 
order corrections in the disorder, as well as the non- 
Gaussian terms neglected in Eq.(^8|) renormalize the ef- 
fective value of ro at small frequencies, so that in Eq. ( |T^ ) 
we should replace ro r{iujm)- Note that r~^{uj) is pro- 
portional to the frequency-dependent conductivity. A 
microscopic calculation of r{uj) is beyond the scope of 
the methods developed in this work. However, we may 
assume different scenarios for r{Lu) and calculate the con- 
sequences for the DOS. 

First of all, suppose that < \imi^^or{uj) = r* < oo, 
so that the system is a normal metal, with a finite static 
conductivity |p^ . In this case the non-Gaussian terms 
neglected in Eq.(p^ are irrelevant (in the renormaliza- 
tion group sense) with respect to the Gaussian fixed point 
action ||T^ . We conclude that if interactions stabilize a 
normal metallic state in 2d |1^, then the DOS exhibits a 
Coulomb gap i^(i^) ~ C|w|/e , just like in the strongly lo- 
calized regime. The numerical value of C depends on r*, 
which we have not calculated. Because the DOS vanishes 
at the Fermi energy, a normal metal in 2d (if it exists) is 
not a conventional Fermi liquid. Let us emphasize that 
for our calculation wc have assumed that the screening 
length is smaller than the mean free path |7j. In 
contrast, in the strongly localized regime the screening 
length is infinite and the Coulomb gap can be explained 
from the classical Hartree energy [Q. Thus, the physics 
responsible for the Coulomb gap is different in both cases. 
Numerical evidence that in 2d the Coulomb gap in the 
DOS survives in a non-Fermi liquid metallic state has 
been found previously by Efros and Pikus |l9| . 

Depending on the low-frequency behavior of r{uj), 
there are two other possibilities. If r(ci;) ^ for w ^ 0, 
the system is a perfect metal. This is the physically most 
plausible scenario, because it is consistent with the gener- 
alized scaling theory Q . In 2d we expect that the diver- 
gence of cr{uj) is logarithmic, so that r{uj) ~ — 7/ln(wTo) 
with 7 > 0. Because the weak logarithmic singular- 
ity does not affect the irrelevance of the non-Gaussian 
terms neglected in Eq. ([l^) , we may calculate i'{uj) by 
substituting ro — 7/ln(wTo) in Eq.(p^). This yields 
i/(w) cx |w|'>'/2. Note that in Ref. d Finkelstein found 
(incorrectly ||l^,|l^) r(cj) ^ — i/ln(cjTo). In this case we 
obtain iy{Lu) cx in agreement with Ref. S. The 
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third possibility is r{uj) oo for a; ^ 0, so that the sys- 
tem is an insulator. Because we know that Eq. is vaUd 
for strongly localized electrons and remains correct 
even if cr(0) is finite, by continuity it is extremely plau- 
sible that i^(ijj) ^ C\uj\/e* in the entire localized regime. 

Experimentally J^(ti^) can be obtained by measuring the 
tunneling conductance dl/dV as function of the applied 
voltage V. Such a measurement was performed in Ref. 
pO| , and the measured dl/dV was compared with the 
perturbative prediction (^. Keeping in mind that a log- 
arithmic dependence was only observed in a rather small 
voltage interval, it seems that the data shown in Fig.l (a) 



of Ref ^ 
Coulom 



2C 



are consistent with the emergence of a linear 
gap (|l|). It would be very interesting to mea- 
sure the tunneling conductance in the 2d materials that 
exhibit a metal- insulator transition j^]. On the metal- 
lic side of the transition, a linear dependence of dl/dV 
on the voltage V would be consistent with a finite value 
of the static conductivity, while a non-linear dependence 
would be consistent with <t(0) = oo. 

I would like to thank S. Chakravarty for many discus- 
sions and for pressing me to sharpen my arguments. I am 
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(c) /\A/ = + r<j~^^^y\/^ 



FIG. 1. (a) This exchange correction to the average 
Green's function gives rise to the In^-correction to the DOS 
in Eq.(^). The solid arrows represent non-interacting disor- 
der averaged Green's functions, the shaded semi-circles denote 
Difi'usion corrections defined in (b), and the thick wavy line 
denotes the dynamically screened interaction defined in (c) 
and (d). Here the dashed line denotes impurity scattering, 
and the thin wavy line is the bare interaction. 
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